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Abstract 



Flow induced by a random acceleration field (g-jitter) is considered in two 
related situations that are of interest for microgravity fluid experiments: the 
random motion of an isolated buoyant particle and coarsening of a solid-liquid 



. mixture. We start by analyzing in detail actual accelerometer data gathered 

during a recent microgravity mission, and obtain the values of the parameters 
defining a previously introduced stochastic model of this acceleration field. 
We then study the motion of a solid particle suspended in an incompressible 
fluid that is subjected to such random accelerations. The displacement of 
the particle is shown to have a diffusive component if the correlation time 
of the stochastic acceleration is finite or zero, and mean squared velocities 
and effective diffusion coefficients are obtained explicitly. Finally, the effect 
of g-jitter on coarsening of a solid-liquid mixture is considered. Corrections 
due to the induced fluid motion are calculated, and estimates are given for 
coarsening of Sn-rich particles in a Sn-Pb eutectic fluid, experiment to be 
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conducted in microgravity in the near future. 



Typeset using REVT^ 



I. INTRODUCTION 



With recent frequent access to a microgravity environment, more attention is being paid 
to a precise characterization of the effective acceleration environment onboard spacecraft 
(g-jitter), as well as to the analysis of potential effects of such an environment on a number 
of experiments, compared to an ideal zero gravity situation j^-^. We study in this paper 
the motion induced on particles that are suspended in an incompressible fluid by an external 
random acceleration field. As an extension, coarsening of a solid-liquid mixture is considered, 
and the effects of g-jitter estimated for the case of a Sn-Pb eutectic. This system will be 
studied in microgravity in the near future. 

Whereas qualitative information on the residual acceleration field onboard spacecraft has 
been available for some time, it is only recently that a systematic effort has been made to 
collect long temporal sequences of acceleration data over a fairly wide frequency range . 
In the frequency range we study (10^^ — 10^ Hz), the SAMS project now routinely determines 
for each mission the three components of the residual acceleration field at selected points 
in the spacecraft. This includes, in some cases, sensor heads at the same location where an 
experiment potentially susceptible to this kind of disturbances is being conducted. 

From a theoretical point of view, the first issue to be addressed concerns the introduction 
of a suitable model of the residual acceleration field. Based on available accelerometer data 
and their associated power spectra, it seems necessary to distinguish between frequency com- 
ponents that can be modeled as systematic (or deterministic), and thus traced back to some 
mechanical device producing a periodic disturbance of known amplitude and frequency, and 
random components arising from a number of independent sources with variable frequencies 
and intensities. We further note that physical sources of accelerations contribute to the 
overall acceleration environment in two ways: directly, and indirectly by exciting some of 
the natural vibration modes of the spacecraft. 

Most of the studies to date have focused on a deterministic acceleration field modeled as 
a superposition of periodic functions of fixed amplitudes and frequencies |7HrTl • Also, some 
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studies have considered the effects of short and isolated pulses 0. The approach presented 
in this paper, on the other hand, models g-jitter as a random process in time |T2|,|13[. We 
assume that the process obeys Gaussian statistics, consistent with the assumption that many 
independent sources contribute to the acceleration field, thus requiring the knowledge of only 
the first two statistical moments of the effective acceleration field g{t). We choose (g) = 0, 
where ( ) denotes an ensemble average. In general, a nonzero average can be incorporated 
into a steady component, and g redefined to be the deviation from the average. We do not 
consider here the effect of this component. To model the statistical behavior of the second 
moment we chose a Gaussian process called narrow band noise, defined by the correlation 
function, 

QAt-t') = {g^{t)g,{t')) = S,,{g')e-\'-'y^cosn{t-t'), (1) 

where gi{t) is any of the three components of the acceleration field, and (g^), t and Q 
are three constants which characterize the process: its intensity, a correlation time and a 
characteristic angular frequency. A particular advantage of this process is that it allows 
interpolation between two well known limits: the white noise limit when Qt with 
{g"^) T = D finite, in which no frequency component is preferred, and monochromatic noise 
when VLt oo, ((7^) finite, in which each realization of the noise is a periodic function of 
angular frequency VL. In this case, the ensemble refers to a distribution of amplitudes and 
phases, with identical angular frequency for each realization. Monochromatic noise is akin 
to the deterministic studies in which g{t) is modeled by a periodic function, but still retains 
random values of the amplitude and phase. 

In Section p, we present a statistical analysis of a long time series gathered by the 
SAMS team during the recent SL-J mission. A window of approximately six hours (sampled 
at 250 Hz) is analyzed to determine the existence of deterministic and random components, 
and to calculate the values of the parameters needed to characterize both. A given time 
series can appear to be deterministic or stochastic depending on the range analyzed: If 
a random function is correlated over times of the order of r, its time series will appear 
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deterministic when analyzed over time windows T <^ t and random otherwise. In the 
case of the SAMS time series that we have analyzed, we find that there is a systematic 
or deterministic component of frequency 17 Hz. The rest of the spectrum is comprised 
of a superposition of random components with small correlation time, and a white noise 
background (with a correlation time no longer than the sampling period of 1/250 s. We 
also find significant deviations from Gaussianity, mainly in the larger amplitude impulses. 
It may be necessary to introduce other stochastic models that are not Gaussian to study 
these contributions (shot noise, for example). 

Section |T| considers the motion of a particle suspended in an incompressible fiuid of 
different density, when the fluid is subjected to the acceleration fleld described in Section 
|I[ If the residual acceleration fleld is deterministic and periodic, the suspended particle 
performs an oscillatory motion, with both velocity and displacement bounded. If, on the 
other hand, the acceleration fleld is random, the mean squared velocity of the particle is 
bounded, but its mean squared displacement grows linearly in time. The effective diffusive 
coefficient is obtained as a function of the parameters of the fluid and noise. In particular, 
we note that measurement of the mean squared displacement of a suspended particle in a 
microgravity environment would provide an independent data set from which one could infer 
the parameters that characterize the residual acceleration fleld. A similar principle has been 
used to design a passive accelerometer system fl^ to obtain the steady component of the 
residual acceleration fleld from the linear drift of a suspended particle. Not surprisingly, it 
has proved difficult in actual microgravity conditions to maintain a well defined alignment 
of the container with respect to the residual acceleration field. 



Section |^ addresses the effects of g-jitter on coarsening of a solid-liquid mixture. A 
random acceleration field induces a random velocity field that may lead to enhanced coales- 
cence and solute transport. The analysis focuses on the solid-liquid mixture Pb-Sn which 
will be studied in microgravity in the near future. We find that g-jitter effects are small for 
the conditions of the experiment, and with the values of the noise found in Section 
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II. TIME SERIES ANALYSIS OF G-JITTER DURING THE SL-J MISSION 



The analysis described in this section is based on actual g-jitter data collected during 
the SL-J mission (SAMS-258), that flew on September 13-20, 1993. We have focused on 
the head A SAMS detector, and studied the series during the time window MET 0017 to 
MET 0023, roughly a period of six hours. All three Cartesian components of the residual 
acceleration field have been included in the analysis. The sampling frequency is 250 Hz. 
The data used was gathered continuously throughout the period mentioned, with automatic 
re-calibration of the sensor heads when needed (corrections for the signal gain have been 
taken into account according to the calibration data also gathered during the mission). 

We do not focus here on some basic statistical properties of the signal which are already 
automatically monitored (its running mean and root mean square values), but address two 
basic points: (a) the existence of deterministic and random components during this partic- 
ular observation period, and (b) the Gaussian nature of the time series. 

Consider first a temporal series g{t) and its power spectrum over a finite window [— T, T] 
defined by 

PT{n) = ^J^^C{t)e-^"-¥'dt, (2) 

where C{t) is the autocorrelation function defined in Eq. (|l]). The autocorrelation function 
can be obtained from Prin) as a Fourier series, 

oo 

CT{t)= E PT(n)e^^, (3) 

n=— oo 

where we have introduced the notation CT{t) to indicate that Crit) = C{t) in [— T, T], and 
is periodic outside of this interval. The integral in Eq. @) can be evaluated exphcitly to 
yield, 

^e-^^ (-A cos(^ - n)T + sin(^ - fi)r) + A 



A2+ -1]]' 



+ 



^ -A cos(f + n)T + (2^ + fi) sin(f + n)T) + A 



A2+ (f + r] 
6 



(4) 



where A = 1/r. In the white noise hmit Prin) = D/T, whereas in the monochromatic hmit 
mT/T,n + ^nn/T-n}- In the first case, the dominant contribution comes from 
the term 1/r in Eq. (|^), whereas in the second case it comes from the term proportional 
to sin(-^ — Q)T. The important point to notice is that in the white noise hmit Prin) is 
inversely proportional to the window size T, whereas in the monochromatic limit, Prin) is 
independent of T. Thus, we argue, an analysis of the power spectrum Pr(n) as a function 
of the window size can provide information on the existence of deterministic or random 
contribution, at least within the available ranges of T and r. 

These results are in fact more generally valid and not restricted to narrow band noise. 
Consider the integral, 

gT{n) = ^f_^dte-^'^g{t). (5) 

If g{t) is a random process, with a correlation time t T, then for each n grin) is the 
sum of approximately 2T/r statistically independent variables. Therefore, according to the 
Central Limit Theorem, the integral will obey Gaussian statistics, with variance 0(T/t). 



As a consequence, grin) ~ 0{j;^T /r) = 0{^) or Pxin) = {gxin)]'^ ~ ^ij^)^ iii agreement 
with the result obtained for narrow band noise in the limit of short correlation time. On the 
other hand, for most deterministic functions grin) ~ (9(1) instead, and Prin) is independent 
of the window size T. 

We have obtained an estimate of Prin) for the time series of g{t) obtained during the 
SL-J mission and for a range of values of T. Since the time series is discrete, we consider 
windows comprising N data points, with NAt = 2T where 1/At = 250s~^ is the sampling 



rate (further details on various methods to estimate power spectra can be found in ref. [0). 
Briefly, the power spectrum for a stationary process, one in which its statistical properties 
are independent of time, is calculated by averaging Pj-, also known as the periodogram. 
The relative statistical error associated with a single periodogram is 100% for all frequen- 
cies. Reduction in error by a factor proportional to 1/Np can be achieved by averaging Np 
periodograms calculated over disjoint time intervals. The estimate of the power spectrum 



presented here is obtained from approximately 6 hours of data, sampled at 250 Hz. Each 
periodogram is calculated for a fixed number of sample points (beginning with 64 and in- 
creasing by factors of two) and then averaged over the entire 6 hour period. The resulting 
estimates of the power spectra are summarized in Fig. |1]. The power spectrum is broad- 
band, with a few peaks at fixed frequencies. The background intensity does decrease with 
increasing N, indicating its random nature. 

To further elucidate the scaling with T, we show in Fig. ^ the value of Pt at selected 
frequencies as a function of T. The frequency components shown in this figure include 
the peaks of Fig. |l|, and one intermediate value. Three types of behavior emerge. First, 
the value of Prif = 17) Hz is independent of T for the range of window lengths studied. 
Therefore, and within this range, this component appears to be deterministic in nature with 
an amplitude ^ < > = 3.56 x lO^^^g^, where qe is the intensity of the gravitational field 
on the Earth's surface. There are two additional components that have a finite correlation 
time. We have fitted the amplitude of the peak to {g'^) ^ (l — e^^/"^^ /2T and estimated 
for the component at 22 Hz y/< g'^ > = 3.06 x lO'^^g^ and r = 1.09s, whereas for AAHz 
we find \J < g'^ > = 5.20 x 10^'^gE and r = 0.91s. As an estimate of the white noise 
background, we obtain from the slope of the intensity of the 8Hz component versus N the 
value D = 8.61 x 10~^cm^/s^. 

In summary, assuming that the various frequency components can be studied indepen- 
dently, and that they are independent of the broadband background, we conclude that the 
time series analyzed contains a deterministic component (i.e., a component with a correlation 
time larger than the largest window studied), a few isolated components of large amplitude 
but small correlation time, and a fairly constant background, of smaller amplitude, and very 
small correlation time. 

To further investigate the statistical nature of the acceleration we calculate a number of 
statistical moments of g{t). We first present the (one-point) probability distribution of g 
obtained from a histogram of the time series. The histogram comprises 200 bins of width 
O.OOSgE- The result is shown in Fig. (||), together with a fit to a Gaussian distribution. It 
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is apparent that the distribution is substantially of Gaussian form at low amplitudes, but 
there are significant deviations near the tails. We only show in the figure the x component of 
the acceleration field. The distribution for the other two components is virtually identical. 

Figure (^) presents the results of higher statistical moments. Normalized cumulants have 
been introduced as follows: 

_ « gjorgitr » .g. 



where << ... >> is the standard cumulant For a Gaussian process, all cumulants 

should be zero except for C2o(t) = Co2(t) and Cu(t). Also note that Ci2(0) = C2i(0) reduces 
to the standard definition of skewness of a distribution, and C22(0) to its kurtosis. Again, 
significant deviations from Gaussianity are found. Further analysis is needed to elucidate 
whether the deviations from Gaussianity in both Figs. and originate entirely from the 
deterministic component at 17 Hz, or are a more intrinsic feature of the random components. 

III. MOTION OF A BUOYANT PARTICLE 

We discuss in this section the motion of a suspended particle in an incompressible fluid 
of different density, when the fluid is subjected to an effective acceleration field of the type 
described in Section This type of motion has also been termed inertial random walk, 
because of the similarity with Brownian motion. The difference, of course, is that the random 
motion of the particle is not induced by thermally induced collisions with the molecules of 
the fluid, but it results from an effective random buoyancy force acting on the particle. A 
qualitative analysis of this process was already given in . 



Consider a spherical particle of radius R and density pp submerged in an incompressible 
fluid of density pf. If the fluid is enclosed by perfectly rigid boundaries, the buoyancy force 
acting on the submerged particle is Ff, = ^tt {pp — pf) R^g(t), where git) is the effective 
acceleration field. In the frame of reference co-moving with the container enclosing the fluid, 
g{t) is a body force, with intensity equal to the value of the acceleration of the container. 
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For containers of reasonable size in a microgravity environment, g can be assumed to be 
spatially uniform. Viscous friction will act on the particle. Neglecting memory terms and 
corrections due to the finite size of the container, the viscous force is given by Stokes' formula 
= —QnriRv, where t] is the shear viscosity of the fluid, and v the velocity of the particle 
relative to the fluid at infinity. For simplicity, we consider in what follows a one dimensional 
case and write, 

x + 'yx = Apg{t), (7) 

with 7 = Qt]/ {2ppE?) and Ap = [pp — Pf)/ Pp- This is a standard Langevin equation which 
can be integrated formally to solve for the velocity v = x. Squaring the formal solution and 
taking the ensemble average one finds, 

•J 

Ap^e-^^* /* dt' f rft"e^(*'+*") {g{t')g{t")) . (8) 
Jo Jo 

The first two terms account for the statistics of the initial condition and are therefore 
arbitrary. Note that in general, the correlation between v{t) and the noise g{t') for t' > t 
does not vanish provided that g{t) has a finite correlation time. For monochromatic noise 
(r — s> oo), this correlation is a non-decaying oscillating function of the time difference {t — t') 
in the limit t,t' —>■ oo. Therefore, the second term of the right hand side decays at least 
as e~^^ and does not contribute for 7^ ^ 1. Two cases of initial conditions are particularly 
simple. First the case when the particle starts at t = with zero velocity and has not been 
subject to the influence of the noise for t < 0. In this case (f ^(0)) = and {v{0)g{t))=0 and 
Eq.(ll) reduces to 

,2\ ^P'{9')r 



V 



(7 + A)' + fi2j ^(^ _ A)^ + f]2 

(1 - ^ (A^ - 72 + n') + (1 + A (-A^ + 7^ + fi^ 

+2A (A^ - 72 - l]2j g-(7+A)t _ 4nA2e-(^+^)* sin nt] . (9) 



\2 



X 
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The second and physically most relevant choice for the initial condition is to take precisely 
the statistics of the steady state. In this case, the transient terms of Eq.(^ will cancel 
exactly the decaying terms of Eq. giving a constant value for (f ^) 

_ Ap2 (^2^(^ + 1) 



oo 



7 



(10) 



(7 + ^) 

This is also the solution for arbitrary initial conditions when 7^ ^ 1. 

Consider now the limits of white and monochromatic noise. For white noise, one has, 

T— »0 \ / 00 'y ^ 



and in the monochromatic case. 



ii-(-L=y^- (12) 



r — ^00 \ / 00 



(An average over phases of the deterministic forces is assumed in the ensemble average of the 
monochromatic noise limit, otherwise {v"^)^ would be an oscillatory quantity.) In the white 
noise limit, {v"^)^ is given by a fluctuation-dissipation relation since Ap^D is the intensity 
of the fluctuations, and 7 the intensity of the dissipation. In the monochromatic limit, 
however, (f^)^ oc 1/7^ (for low frequencies, Q/'y <^ 1). This is precisely the overdamped 
limit of Eq. (|^). In all cases the mean squared value of the velocity saturates at a finite 
value at long times. In the monochromatic noise limit, this is the case even in the limit of 
small viscosity 7 — > 0. In the white noise limit, on the other hand, viscosity is essential for 
saturation. 

The equation for (x(t)^) can be obtained by multiplying Eq. ([^ by x, and rewriting it 

as. 

This equation is cumbersome to solve in the general case of narrow band noise. We focus 
separately on the two limits of white and monochromatic noise. In the standard case of 
white noise, we have {x{t)g{t)) = 0, and {v{0)g{t > 0)) = and (from Eq. (|])), 
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Ap^D 

7 



(l - + ^^(0)2) 



(14) 



If we assume (f (0)^) to take the steady state value ^^'^ (or equivalently, for an arbitrary 
initial condition if 'jt ^ 1), and taking the ensemble average of Eq. ([I3l), we have for 7t ^ 1 
{x'^{t)) = ^^^a ^ t + constant, that is diffusive motion with an effective diffusion coefficient 

For the monochromatic case, the displacement of a particle with x(0) = and v{0) = 
can be easily obtained. Square the deterministic solution for g(t) = gocos{Qt + (p) and 
average over the phase 0. With the identification {g"^) = gl/2 and for 7^ ^ 1, we obtain 



(72 + fi2)2 



1 



7 



+ 



7 



7 



cos Vtt + sin fit 



• (15) 



In this case the mean square displacement of the particle is bounded. For the particular 
case of 7 = 0, Eq. (^) is not valid and the mean square displacement of the particle may 
be unbounded. 

It is worth remarking that Eq. ( |l5l) is only valid in the strict limit r = 00. If the corre- 
lation time T is finite, there is always a diffusive behavior superimposed to the oscillations. 
An effective diffusion coefficient -De// can be defined as. 



1 /•*+^ 1 
lim — 



DeSf 

For narrow band noise we find, 
lim {x{t)g{t)) = Ap{g^ 
and hence 

Deffl 
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v\t')) + Ap{xit')git')) dt', 



T = — . 
fi 



(16) 



A(A + 7) -fi2 



((A + 7)2 + fi2)(A2 + fi2) 



^(A2-72 + fi2)+^2_p^^2 



+ 



+ Oscillatory terms. 



A(A + 7) -fi2 



(17) 



(18) 



Ap2((^2) ((^ + A)2 + fi2)((^_ A)2 + f]2) ((^ + A)2 + fi2)(A2 + (]2)- 

As a first correction to the monochromatic limit, the effective diffusion coefficient to first 
order in 1/r reads 



I? 



eff 



Ap' {g') 

(72 + Q^y 
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2 + -T + 7^ - + 0(l/r2). 



fi2 / r 



(19) 
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IV. COARSENING IN SOLID-LIQUID MIXTURES 



We consider next coarsening of a solid-liquid mixture being subjected to a fluctuating 
acceleration fleld of the type described in section |T|. Such a study is relevant in connection 
with an experiment that will be conducted in microgravity in the near future. The exper- 
iment will address the asjTiiptotic power law growth governing coarsening, as well as the 
dependence of the amplitude of the growth law on the volume fraction of the precipitate 
phase. The absence of gravitationally induced sedimentation will allow a careful quantitative 
study of these two important theoretical issues. 

A residual acceleration fleld can produce a number of deleterious effects on otherwise 
purely diffusive controlled coarsening, which we address in this section. We focus here on two 
such effects: random motion of the suspended particles induced by the effective (random) 
buoyant force and the concomitant increase in the likelihood of particle coalescence, and 
additional flow in the fluid phase caused by g-jitter and its effect on solute mass transport. 

In order to obtain numerical estimates for these two effects, we will consider experimen- 
tal parameters for a solid-liquid mixture of Sn-rich particles in a Pb-Sn eutectic liquid, the 
system that will be used in the microgravity experiment ||T8|JT9|1 . The density of the pre- 
cipitating solid phase is Pp = Y.OSSg/cm^, whereas that of the liquid is pi = S.OTAg/cm^. 
The kinematic viscosity of the liquid is u = 2.48 x 10^'^ cm^/s, and the solute diffusivity is 
Ds = 5.6 X 10~^cm^/s. It is anticipated that coarsening will be studied for a period of 5 
hours, with an average particle size at the end of that period of Rav — 7 x lO^^cm = 70 pm. 
Given the size of the particles and the small values of the residual gravitational fleld, in- 
ertial effects will be completely negligible. In what follows, we focus almost exclusively on 
Stokesian dynamics for the suspended particles pO[. In addition, the solution is not mono- 
disperse, but rather a scale invariant particle size distribution evolves dynamically , with 
particles sizes ranging from to ~ 1.5Rav In all the estimates that follow, the average 
size at the end of the flve hours is used, an overestimate for most of the duration of the 
experiment, and a slight underestimate at the latest times. As will be seen below, a factor 
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of two in Rav would not modify our conclusions. 

Neglecting inter-particle hydro dynamic interactions, precipitate particles will execute a 
random motion of the type described in Section For the case of monochromatic noise 
(fixed frequency and random phase), the average quadratic displacement of each particle 
remains bounded, and is given by Eq. (|T3]). For the parameters of the fluid given 7 = 260s~^, 
and by using the amplitude of the 17 Hz component of the power spectrum in Section 
H (^] = 27r 17s-\< / >= 1.27 x 10"^^|), we find that max {< >} ^ IQ-^cm^, 
and hence negligible. At the other extreme, we find that for white noise (replacing the 
factor of 2 by 6 corresponding to diffusion in three dimensional space in the expression 
for the effective diffusion coefficient), the mean squared displacement after five hours is 



< > (t = 5 hr.) = 8.85 x lO^^cm? or v<x^> ~ 30/im. Clearly the average square 

displacement induced by the white noise component of the residual acceleration field is 

much larger than that induced by the monochromatic component, but it is still about one 

half of the average particle size. Therefore random motion of particles induced by g-jitter 

will not lead to significant coalescence during this time period. As a reference, we quote the 

average squared displacement induced by thermal Brownian motion: (x^) = kBTt/3'7iRavV 

or (x^) ~ 39 /im, for the same time span and taking Rav = 70 /im. Both effects are 

therefore expected to be of the same order of magnitude during the experiment. 

It is also possible to estimate hydrodynamic interaction effects between pairs of spherical 

particles to show that it leads to weak attraction at long distances, and repulsion at short 

distances. The relative displacement r of particle 2 with respect to particle 1 immersed in 

an incompressible fluid satisfies p^j23[] , 

dv^ — * — * 

— = (cjai - i^n) ■ Fi + (0^22 - ^^21) ■ F2, (20) 
where Fi is the force acting on the i-th particle, and ojij are hydrodynamic mobility tensors. 



given, e.g., in references [|22| , |23| . After some straightforward algebra, the leading contribution 



at distances large compared to the particle radii is given by 
df _ 2{pp - Pf) . 2 o2\ , (Pp ~ Pf) - ^2) 1 



{Rl - Rl) m 



dt 9/1 ^ ' 3/i 
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m, (21) 



where X is the identity tensor. The first term in the right hand side describes the relative 
motion of two independent particles of different size, and therefore its magnitude has already 
been estimated above. Both the longitudinal and transverse components of the second term 
in the right hand side of Eq. (|2T| ) are of the form, 

where, for the longitudinal component, A= {pp — Pf){R\ — Rl)/3fi. 

Consider an initial inter-particle separation Ri. In this case, and for times shorter 
than the average time needed for the two particles to coalesce, the quantity y = jIA is a 
Wiener process if g(i) is Gaussian and white, and therefore the conditional probability for 
r is, 

P(r,t|ro,to) = ^ e isj^a^c-'q). (23) 

|A|^47rD(t-to) 

The ensemble average of r, < r > can be computed analytically, 



< ^ >= V ^ [2^(^ - ^o)] e"32^^^(*-o)D_3/2 ( , 1 , (24) 

V 4 I \ 2\A\J2D{t-t,y 



where Dp{z) is a parabolic cylinder function (formula 9.240). For short times, the 
asymptotic form of Dp{z) for large z allows the computation of < r > 

/ A'^D(t-to)\ , , 

<r>=r,\\ , (25) 

which decreases in time regardless of the sign of A. Therefore g-jitter induces an effective 
hydrodynamic attraction between pairs of particles. However, for the experimental values 
given above, and taking Rx = 1.5Rav and R2 = 0.5Rav, A = 1.8 x 10~^cm s. If tq ^ 200yum, 
then after 5 hours the inter-particle separation would have decreased by approximately 7/im, 
and therefore small compared to particle radii. 

The attractive interaction is not confined to short times, but it arises directly from the 
1/r dependence in Eq. (P^. By taking the average of Eq. (P^, using the Furutsu-Novikov 
theorem [^] and the fact that the noise is Gaussian and white, one finds. 
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^^^Aol'-im, (26) 

where S/Sg{t) stands for functional derivative with respect to g. Directly from Eq. (0), we 
find that 6{l/r{t))/6g{t) = —A/r^, and therefore, 

d <r > ,o_/l 



dt -^'"{7')' 

identical to Eq. (|25|) with 1/rQ replaced by < >. Since r is a positive quantity, 

d < r > /(it < for all values of r. It is also interesting to note that the effective attractive 
interaction is not confined to the term proportional to 1/r in the hydrodynamic mobility, but 
that attractive contributions arise from higher powers of 1/r as well. In fact, this attraction 
is generic for overdamped motion and multiplicative noise provided that the mobility is a 
decaying function of the inter-particle separation ||26[ . 



The question naturally arises as to the behavior of pairs of particles near contact, or of 
particles near a solid wall. In either case, lubrication theory allows the calculation of the 
mobility tensor. The longitudinal component vanishes linearly with inter-particle distance 
whereas the transverse component becomes non-analytic (diverges logarithmically at short 



distances) |]23|. In both cases, the mobility increases with inter-particle separation leading 
to an average repulsion {d < r > /dt > 0) following the same arguments given above. 

Estimating the effect of g-jitter on mass transport in the fluid phase and therefore on 
coarsening kinetics is far more complex, and we will not attempt a complete solution here. 
We show below that the order of magnitude of this contribution to coarsening is also small, 
and therefore a detailed calculation is not necessary. However, and in order to motivate 



the analysis that follows, let us define an effective Peclet number as Pe = y<v?>R/ Dg, 
where u is the characteristic velocity of the fluid. (The velocity u is of the order of the 
velocity of the particles because the motion of two-phase interfaces due to phase change is 
small in the time scale of change of g{t)). We find for monochromatic noise (Eq. (|T2|)) that 
< >= 3.48 X IQ-^cm^/s^ or Pe = 0.23. For white noise, on the other hand, according 
to Eq. (PT|, < >= 6.39 x 10"^cm^/s^, or Pe = 0.32. Therefore it would appear that 
convective transport of mass is not negligible in front of diffusive transport. 
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Such a calculation, however, overestimates convective transport. Since the system is 
statistically uniform and the effective acceleration field averages to zero, the average velocity 
of the fluid has to be zero. Hence no overall convective motion would result in the (longer) 
time scale over which diffusive transport occurs, once the fluctuating component is averaged 
over times much larger than the correlation time of the noise. Since < m >= but < 

>7^ 0, the lowest order contribution to transport due to the motion of a fluid element 
is diffusive [^]. The calculation of this effective diffusivity can be carried out in a mean 
field approximation. Consider a single solid particle immersed in the fluid phase, so that the 
presence of the remaining particles can be subsumed in a far field composition Coo , and a 



far field velocity Moo? both to be determined self-consistently for a given particle distribution. 
The cut-off distance is typically of the order of the inter-particle separation. 

Fluctuations in Uoo due to the motion of the ensemble of particles (again, fast compared 
to coarsening times) lead to an effective increase in diffusive mass transport, and hence to 
an increased diffusivity. For a quiescent system, the average velocity of the fluid Uav, and 
the average velocity of the particles Vav are related by, 

(t)Vav + (1 - (fyUav = 0, (28) 

where is the volume fraction of the system, which we assume to be small. The subindex 
av indicates an average over the particle distribution at fixed time. In mean field, we take 
Moo = Uav, and also 

-^^2<M^zi;£M + o(^), (29) 
with hydrodynamic interactions contributing to 0{(f)). Therefore, 

u^ = — + C»(0). (30) 



The far field effective diffusivity is then [28 



D,jf = Ds + - (Moo(t) ■ Uooit + 1')) dt' = Ds+ ^ I (31) 
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in the white noise hmit. For the parameters of the experiment given above 

—L— ^ 1 

and therefore numerically negligible. 

In summary, even though anticipated Peclet numbers based on the scale of the flow are 
of order unity, mass transport due to convection is expected to be negligible during the 
solid-liquid coarsening experiment. Since the time scale of acceleration variations is short 
compared to coarsening times, and leads to zero average velocity, the contribution from 
g-jitter to mass transport is diffusive and leads to a very small correction to the solute 
diffusivity. We note, however, that the correction is proportional to l/7a„ oc i?^^ and hence 
it increases quickly with the average particle size of the precipitate phase. Therefore, either 
under different experimental conditions, or in the strict asymptotic limit of very long times 
(and hence large Rav), transport due to transient accelerations would dominate molecular 
diffusion leading to a different asymptotic growth law for the average particle size. 
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FIGURES 
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FIG. 1. Power spectrum as a function of frequency averaged over a six hour interval during the 
SL-J mission. The curves shown correspond to spectra calculated over windows of size (from top 
to bottom) N = 64, 512, 4096 and 32768. The amplitude of the peak at / = 17 Hz is independent 
of N, whereas the amplitude of the peaks at both 22 Hz and 44 Hz decreases with N. Also in 
this latter case, the shape of the peaks is independent of N. 
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FIG. 2. Amplitude versus window size A'^ for a few selected frequencies to display their deter- 
ministic or random nature. The amplitude of the 17 Hz component remains independent of N 
indicating its deterministic character for the range of window sizes analyzed. Two other compo- 
nents display mixed behavior, with a finite correlation time of the order of 1 s. There is also a 
clear white noise background, exemplified by the amplitude of the power spectrum at 8 Hz. 
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FIG. 3. Histogram of amplitudes of the residual acceleration along a particular direction, and 
a fit to a Gaussian distribution. The distribution is nearly Gaussian for small values of but 
deviates significantly near the wings. Values of the gravitational field intensity are given relative 

to QE- 
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